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SUMMARY 


A review of some approximate modeling techniques used to develop 
analytical solutions for damaged, fiber-reinforced composite materials 
is presented along with previously unpublished results of the past 
year's research concerning the application of these methods to particular 
problems. 

The classical shear-lag stress-displacement assumption is funda- 
mental to much of this work. Based on this approximation, solutions 
are developed for the two-dimensional region containing unidirectional 
fibers with initial damage in the form of, respectively, a notch, a 
rectangular cut-out, and a circular hole. An ultimate stress failure 
criterion is used for both the fibers and the matrix; simple tension 
for the fibers and shear failure for the matrix. Models which account 
for longitudinal matrix yielding and splitting as well as transverse 
matrix yielding and fiber breakage as a function of initial damage, 
material properties and applied stress are presented. The fibers are 
taken as linearly elastic and the matrix material as either elastic- 
perfectly plastic or elastic-strain hardening. A cover sheet con- 
straining the unidirectional laminate is also introduced although 
only its influence on the unidirectional laminate is modeled; the 
stresses within the constraint layer are not computed. 

For ductile matrix composites (boron/aluminum) the results 
indicate that longitudinal matrix yielding and transverse notch 
extension are the most significant forms of damage to Include in order 
for the model to agree with experimental results. The extent of the 



stable transverse damage is shown to be approximately constant, 
independent of initial notch length. Including a cover sheet and/or 
strain-hardening matrix has minor influence. In the case of brittle 
matrix composites (graphite/epoxy) longitudinal splitting is shown 
to be the dominant form of damage. 

Very little difference is found between the results for the 
three types of initial damage, i.e.,the notch, rectangular cut-out, 
and circular hole. In all cases, the presence of additional damage 
changes the nature of the stress distribution in the unbroken fibers. 
For the original Hedgepeth problem of a notched laminate the stresses 
decay as the square root of the distance from the notch tip; including 
longitudinal or transverse damage significantly reduces the stress 
concentration and gives a much more uniform stress state in the un- 
broken fibers. It is shown that this behavior cannot be accounted 
for by introducing an effective notch length or crack tip damage zone 
with a square-root behavior. 

The formulation of the problem for an edge notch in a uni- 
directional half-plane with no additional damage is also developed 
and the appropriate equations are recorded but numerical results are 
not given. This solution forms the basis for the general problem of 
adjoining half-planes of different fiber-matrix properties now being 
developed by the writers. 
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INTRODUCTION 


A major portion of the writer's research over the past few years 
has concerned the development of suitable analytical techniques for 
predicting the stress state and fracture behavior of damaged composite 
laminates. Much of this work has dealt with approximate solutions 
based on discrete fiber-matrix material models with particular simpli- 
fying assumptions used to relate fiber and matrix stresses to fiber 
displacements. The resulting solutions are not complete solutions 
to the equations of elasticity for the two phase region and a sig- 
nificant portion of the study has been the investigation of the 
agreement of the results with experimental data. 

In past reports and technical papers concerning this work, 
specific models and results have been presented in each paper, but 
no unified report has been written. Rather than add one more paper 
covering only the work of the past year it was felt to be an appro- 
priate time in the development of these methods to review the basic 
assumptions and discuss the significance of the models and the results, 
both for the convenience of having a more complete record and to be 
able to make some important observations concerning the nature of the 
different models, about which we have only recently become aware. 

The initial work in modeling a unidirectional composite contain- 
ing broken fibers was presented by Hedgepeth in [1] for broken fibers 
with no longitudinal or transverse damage other than the initial notch. 
This work was extended by Hedgepeth and Van Dyke in [2] for the special 
case of one broken fiber with matrix yielding parallel to the fiber and 


in [3] for one broken fiber with longitudinal splitting in the matrix. 

In all these studies the fiber breaks were assumed to lie on a transverse 
line and the shear-lag assumption was used. One very important feature 
of the shear-lag model is that it simplifies the equilibrium equations 
by removing the transverse displacement dependence from the longitudinal 
equation, and the fiber stress and the matrix shear stress can be 
determined without solving the transverse equation. In reference [2] 
Hedgepeth and Van Dyke used the same model to develop the solution for 
broken fibers in a three-dimensional unidirectional composite con- 
taining broken fibers with no additional damage. A detailed discussion 
of similar material modeling techniques and some simplified solutions 
are presented by Zweben in [4] and [5]. Eringen and Kim in [6] developed 
a modified solution for the original Hedgepeth problem, [1], by extend- 
ing the shear stress-displacement relation to include fiber bending as 
well as axial displacements; no additional damage was accounted for in 
[6]. The stress-displacement relation assumed in [6] allows for better 
satisfaction of a stress free crack surface, i.e. the shear-lag 
assumption does not have sufficient freedom to remove the shear 
stresses from the crack surface. However, Eringen 's model does couple 
the axial and transverse equilibrium equations and gives a somewhat 
more complicated set of differential -difference equations. The 
inclusion of matrix damage and transverse fiber breaks (notch extension) 
in Eringen's model appears to be much more difficult than in Hedgepeth's 
model and such modifications have not been attempted at this time. 

The initial damage in all of the above studies consisted of 
broken fibers in the form of a notch (crack). Franklin in [7j and 
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Kulkarni et.al. in [8] investigated the case of a circular cut-out 
containing no additional damage and found that the stress concen- 
trations in the unbroken fibers was changed very little over the 
corresponding solution for a notch. 

Goree and Gross in [9] extended the Hedgepeth solutions to 
include longitudinal matrix yielding and splitting for an arbitrary 
number of broken fibers and in [10] developed a solution using the 
Eringen model of [6] for a three-dimensional unidirectional composite 
containing broken fibers but without matrix damage. The results of 
[9] gave very good agreement with experimental results for brittle 
matrix composite and reasonably good agreement for ductile matrix 
composites. The inclusion of transverse stable notch extension to 
this model is shown below to make a very significant improvement in 
the ability of the model to represent the behavior of a ductile 
matrix (boron/aluminum) laminate. 

Over the past year solutions have been developed for the following 
problems: 

1. transverse notch extension 

2. constraint (or cover) layer 

3. strain-hardening matrix 

4. rectangular (rather than slit) initial damage region 

5. circular initial damage region 

6. formulation of the edge crack problem. 

In the first section presented below the formulation of the original 
Hedgepeth problem [1] will be presented, using Fourier transform methods 
rather than the influence function technique as used by Hedgepeth as 



this is the foundation for the remaining solutions. The next sections 
will then consider the development of the equations for, respectively, 
longitudinal splitting and yielding of the matrix, transverse notch 
extension, constraint layer, strain-hardening matrix, rectangular 
damage region, circular damage region, and the edge crack solution. 
Results and comparisons of the various models will then be presented. 

Use of trade names or names of manufacturers in this report does 
not constitute an official endorsement of such products or manufac- 
turers, either expressed or implied, by the National Aeronautics and 
Space Administration. 
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FORMULATION OF MATHEMATICAL MODELS 


A . Two - p i mensional Shear-Lag Model with Broken Fibers 

Consider a two-dimensional unidirectional lamina containing 
broken fibers as shown in Figure (1). The development presented in 
this section will be for no additional damage other than the broken 
fibers. This solution was first presented by Hedgepeth in [1]; however, 
as it is the fundamental solution on which all of the following models 
are based, it will be included in this report for completeness. Further, 
the method of solution appropriate for the various extensions to this 
basic model is somewhat different than that of Hedgepeth in [1], i.e., 
Fourier transform techniques are used directly in the present work 
while Hedgepeth devel oped< the solution by means of influence functions. 
For no damage other than the initial notch the two methods are equiva- 
lent. However, for the extension to matrix yielding and splitting at the 
end of the notch containing more than one broken fiber, the Fourier 
transform method is more direct. 

The formulation given below will also develop the solution for 
the matrix normal stress which is not given in [1]. It should be 
noted that it is often attributed to the shear-lag solution that only 
shear stresses exist in the matrix. This need not be imposed on the 
solution, although a fundamental property of the shear-lag assumptions 
is that the differential -difference equations for axial and transverse 
fiber displacements uncouple such that the axial displacement can be 
found without solving the transverse equation. As the fiber stress 
and matrix shear stress are functions of the axial fiber displacement 
alone, these stresses can then be determined without knowing the 


transverse displacements. Once the axial fiber displacement is known it 
is a simple matter to solve the transverse equation for the transverse 
fiber displacement and compute the matrix normal stress between fibers. 

Following Hedgepeth [1], the laminate is modeled as a two- 
dimensional region, shown in Figure (1), having a single row of parallel, 
identical, equally spaced fibers, separated by matrix. The damage is 
taken to consist of an arbitrary number of broken fibers such that all 
breaks lie along the x-axis, but they need not form a continuous break 
(notch). Later in this section the equations corresponding to a notch 
will be developed. It is this solution for a notch that is extended 
in the remaining sections. In all cases it is assumed that the fibers 
have a sufficiently higher elastic modulus in the axial direction than 
the matrix such that the fibers support all the axial stress in the 
laminate. The matrix supports transverse normal stresses and shear 
stresses. 

Admittedly, most unidirectional composites consist of more than 
one lamina with all fibers in each lamina surely not perfectly aligned 
either through the thickness, or within each layer. These variations 
can have a considerable influence on the stress state. For example in 
[11] and [12] it is shown that the shear stress becomes larger as the 
fiber spacing decreases. Local failures may well occur at the critical 
points through the thickness in advance of laminate splitting which 
could give an apparent shear stiffness considerably different from that 
for the matrix alone. It is assumed that such variations can be 
accounted for by an appropriate choice of a matrix shear modulus Gj^^j 
and a shear transfer distance h. It is with this in mind that the 
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following development will be concerned with an equivalent lamina 


where 6j^ and h are to be determined experimental ly for any particular 
laminate. The following fundamental assumptions regarding the stress- 
displacement relations are made: 

1. The axial fiber stress, ^ » in fiber n is given by 

I _ c n 

°F|n '"F dy ’ 

where Ep is the Young's modulus of the fiber and v^ is the axial 
displacement of fiber n. 

2. The shear stress in the matrix, t1 , , between fibers n and 

ln+1 

n+1 is given by 

Gm 

P+1 “IT ^'^n+l “ ''n^ ’ 

where Gj^/h is the equivalent shear stiffness of the matrix and 
v^ is the axial displacement of fiber n. This relation is the 
basic shear-lag assumption. 

3. The transverse matrix stress, a..| , between fibers n and n+1 

I n+1 

is given by 

"Mln+l h ^ n+l"‘'n^ ' 

where Ej^/h is the equivalent matrix transverse stiffness and u^ 
is the transverse displacement of fiber n. 

4. Consistent with the above assumptions it then follows that the 
stress state on a transverse plane is constant between fibers. 
With these definitions and assumptions, the equilibrium equations 


for a typical element as indicated in Figure (2) are as follows: 


Ap dcjp 

4 - T 

t dy n+1 



(A.l) 
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and 


M n+i 



0 . 


(A. 2) 


Substituting from the above stress-displacement relations the 
following pair of differential -difference equations are obtained: 


and 


EpAph d^Vj^ 






+ V.n-2V+V 1 
n+1 n n-1 


= 0 




- 2u + 
n 


'n-1 


+ 


^ A / 

h dy |''n+l 



(A.3) 


(A.4) 


As mentioned above, the equilibrium equation in the axial direction, 
equation (A.3), is seen to be independent of the transverse displace- 
ment and can be solved without solving equation (A.4). Hedgepeth [1] 
does not consider this transverse equation but develops a solution to 
the axial equation (A.3) only. We now proceed to develop a solution 
to equation (A.3) and, with this solution determined it will then give 
the solution to equation (A.4) for the transverse fiber displacements. 

Noting the coefficient of the first term in equation (A.3), the 
following changes in the variables are suggested: 


y = 


A| 


EpAph 




Or- = a 

F n 


- r ^ 

n dy 


(A. 5) 


and 


v = cr 
n ° 


1 


Aph 


EpGpit 


n 


where n , , V are nondimensional . Equation (A.3) then becomes 

h n ^ 

independent of all material properties as 


d^V, 




= 0 


(A. 6) 
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This differential -difference equation can be reduced to a differential 
equation by defining a new function V(n.0) such that the normalized 
displacement V^(n) is the Fourier coefficient in a Fourier series ex- 
pansion. That is, 

00 

V(ti,e) = z V (n)e’''"® , and (A.7) 

n= - oo 

as the displacements are continuous functions of ns this representation 
is necessarily valid and can be inverted to yield 

1 'H' _ ■} n o 

Vn(n) - 2 ^ ^ V(ns6)e de . (A. 8) 

-7T 

Substitution for V^(n) in terms of V(ns0) in equation (A. 6) re- 
sults in 

2[1 -cos(e)]) e'"® de = 0 . (A. 9) 

-7T I dn 

This equation is of the form 

^ / F(ns6)e^'^® de = 0 , for all n and n. 

~TT 

The function F(n,0) is continuous in 6 and therefore, if the integral 
is to vanish for all n the function F(nsO) must be zero. The equation 
specifying V(ns0) is then 
2 - 

^ = 0 . (A. 10) 

dn 

where 

6^ = 2[1 -cos(e)] = 4 sin^(e/2) . 

The solution to the problem of vanishing stresses and displace- 
ments at infinity and uniform compression on the ends of the broken 
fibers will now be sought. The complete solution is obtained by adding 



the trivial results corresponding to uniform axial stress and no broken 
fibers to the following solution. 

The boundary conditions are, 

Vn(n) =0 as n “ j (A. 11) 

for all fibers, 
dV (n) 

5p|^ = - 1 . at T, = 0 , (A.12) 

for broken fibers, and 

V„(n) = 0 , at n = 0 , (A. 13) 

for unbroken fibers. 

Equation (A. 10) has the complete solution satisfying vanishing 
stresses and displacements remote from the damage as 

V(a,e) = A(e)e‘^'' , (A. 14) 

where the function A{e) must be determined from the remaining boundary 
conditions. 

Using equations (A. 8) and (A. 14) the displacement is given by 

Vn(n) = ^/ A{e)e'^V"® de. (A.15) 

-TT 

Noting the form of the above integral for n = 0, the boundary condi- 
tion of no displacement for the unbroken fibers, equation (A. 13) is 
identically satisfied by taking 

A(8) = E B.e’''*'® , (A. 16) 

z 

with z being the index of each broken fiber and the are constants. 

One then has precisely the number of constants as broken fibers. 

These constants B^^ are determined from the remaining boundary condition, 
equation (A.12), of a unit compressive stress on the broken fibers at 
1 0 



n = 0. Substituting equation (A. 16) into (A. 15), equation (A. 12) then 
gives, at n = 0 , 


^ / I da = 1 ; fi = 2 sin(e/2) (A. 17) 

“IT 

for z and n corresponding to all broken fiber indices. This equation 
then gives a system of linear algebraic equations for the unknowns 
and the solution is complete. If the broken fibers are symmetric about 
the zero fiber, and therefore form a notch, the above equations are con- 
siderably simplified. Let the center fiber be located by the index n = 0, 
then equation (A. 7) may be written as a cosine series as 


V (n) “ 

v(n,e) = + Z V„(n)cos(na) , (A. 18) 

n=l " 

from which 


V.(ti) = - / V(ri,9)cos(n6)de , (A. 19) 

IT 0 

where, as in equation (A. 14), V(nse) = A(e)e~'^’^ . 

Equation (A. 16) is a cosine series for this symmetric case and is given 

by N 

A(e) = E B cos(jl0) , (A. 20) 

£=0 ^ 

where N is the index of the last broken fiber. One then has N + 1 un- 
known constants B^ with the solution given by satisfying the boundary 
condition of uniform stress on the broken fibers at n = 0, equation 
(A. 12). This gives 

« N TY 

— E B. / 6 cos(ji0)cos(ne)d0 = 1 , n = 0,1,. ..,N (A. 21) 

£=0 ^0 

where, as before, recall that 6 = 2 sin(0/2). For example, if N=0 


which corresponds to one broken fiber, equation (A. 21) gives directly 
as 

B. = ir/[2 / 2 sin(0/2)de] = f = A( 6 ) . (A. 22) 

U 0 O 

The maximum fiber stress in the first unbroken fiber is at ti = 0 
and is given by 


Op ,(0) dV,(0) . ^ , 

— J = ^ A(e)cos(e)d 6 = 3 . 


or, for a unit stress at infinity and an unloaded free end of the broken 
f i ber 


1 


= 1 + 


'Fll 


( 0 ) 


= 1 + 


i _ 4 
3 3 ■ 


The normalized crack opening displacement is given by equation (A. 19) as 
2 Vq( 0 ) = f . 

For the general case of 2N + 1 broken fibers, and with the axial 
fiber displacement now known, the transverse fiber displacement Uj^(n) 
may be determined from equation (A. 4). For the specific case of a 
symmetric number of broken fibers in which the axial displacement is 
given by an even valued transform, equation (A. 19), the transverse 
displacement will be odd valued and is given by 


where 


U(ris 6 ) = 2 U (n)sin(ne), and 

n=l ^ 


U{ii, 0 )sin(ne)de 

n 7T Q 


Un(n) = 


u„(n) 


1 


Aph 


(A. 23) 


(A. 24) 
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Substituting this representation into the normalized form of equation 
(A. 4), the resulting transverse displacement is then given by 


1 

U„(n) = - - ^ 

n' ' ir 


M ^ 


/ dV(n,9) 

ApEph ^ dn 


sin( 0 )sin(ne) . 

1-COS{0) 


(A. 25) 


In the following sections the above solution for a symmetric notch 
will be modified to account for the various damage models. 



B. Two-Dimensional Shear~Lag Model with Broken Fibers and 

Longitudinal Matrix Splitting and Yielding 

The solution developed in section A will now be extended to in- 
clude longitudinal splitting and yielding of the matrix as shown in 
Figure (3). This solution along with extensive results is given in [9]. 
All the previous assumptions are assumed valid and it is only necessary 
to account for the additional damage region parallel to the fibers. It 
is assumed that splitting and yielding of the matrix initiates at the 
notch tip and progresses longitudinally between the last broken fiber 
and the first unbroken fiber as shown in Figure (3). The matrix ma- 
terial is assumed to be elastic- perfectly plastic. 

The same free-body diagram of section A, shown in Figure (2) is 
considered with additional conditions for the last broken fiber, denoted 
by n = N, and y £ L, to account for the longitudinal damage taken as 

t| = -T< y-i > (B.l ) 

'N+1 

where 

<y-£ > = 1 , y 1 ^ 

<y~i> = 0 , y < £ , and (B.2) 

L equals the total damaged length, £ the split length, and the matrix 

yield stress. Yielding is assumed to occur when the matrix shear strain 

reaches the yield strain, Yq- Splitting occurs at a multiple of as 
given by the particular matrix material. 

The equilibrium equations in the longitudinal and transverse 
directions, respectively, for all fibers n, with the exception of N and 



t dy 


+ T 


„ , 1 - T =0 
n+l n 


(B.3) 


and 


M n+l 


"’Hln + {-^In+l+^ln) “ 


(B.4) 


For fiber N,y £ L, x = -T^<y-^-> , and the equilibrium equations are 


Arr da 


and 




- X <y-a> 
dy 0 


I -4. h d 
N+l “ N 2 dy 


T f, = 0 


Tg<y-s.>+ T 


= 0 


(B.5) 


(B.6) 


For fiber N + l, y £ L, = -TQ<y-J2.>, and the equilibrium equations 

are 


^F ‘^'^F N+l 


dy 


N+2 *^0 ^ ^ ^ ’ 


and 


h d 

"mIn+Z ■ °m|n+ 1 2 dy i'^L+2 ■ ^ ' 


(B.7) 


(B.8) 


Substituting the stress-displacement relations into the equilibrium 
equations, the following pairs of equations are obtained: 

For all fibers except N and N + l when y £ L , 

EpArh d^v. 

and 


pAph d"^v f 


(B.9) 


h” (“n+l ■ ^“n “n-l) 2 ” dy (''n+l ' ''n-1 ) " ° ' (B.10) 


For fiber N,y < L , 


EpAph d v,M . 

— p— z 0 - -I - - j!— X <y-A > = 0 , 

dy2 N G„ o J' 
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(B.ll) 



and 


h” {“n+ 1 ■ ? dy { h” ^^''n " ''o ‘ 


For fiber N + 1 , y < L . 


Er-Ar-h d v»|,T hx^ 

T~ ^^N+2 " ''n+I > 


V 


dy' 




(B.13) 


and 


h” (“N+2 ■ ^“N+1 '^n] 2 dy ( h” ^^''n+2 ''n+1^"'o < > ) -0- (B.14) 


Again, the axial equilibrium equation is independent of the trans- 
verse displacement, u^ . 

The three equilibrium equations in the longitudinal direction are 


then: 

for all fibers, except N and N+1 when y £ L , 


EpAph dS 

“G^r -TT* VT^Vn + v^., = 0 , 
M dy 

for fiber N, y < L 


(B.15) 


^cAph d V|^ . 

— p - - 1 — — 7^ Vm T - Vfj - 77— T_ < y-5, > = 0 , 

V dy2 N G„ 0 


(B.16) 


and for fiber N+1, y < L 


EpAph d 

J^2 ''N+2 ■ ''n+ 1 Gj;^ ■''o " 




(B.17) 


The same change in variables as before will be used, with the following 
additional terms: 


T = G 

0 " 


Epht "^0 


1 6 


(B.lSa) 



(B.18b) 


dV. 


Oc ~ o 
F n 


dn 


^0 


Epht dV^ 


Vf 


dn 


T = a 

n 


(»n-Vl)->{«„-V,} 


(B.18c) 


L 



EpAph 


and 



EpAph 


(B.18d) 


where n , a , V , , a, and 3 are nondimensional . 

n n ’ 0 

In these equations, Ep , Ap , t, L, and z are taken as actual fiber 
modulus, fiber cross-sectional area, lamina thickness, and damage dimen- 
sions, respectively. The quantities Tq and Gj^/h are equivalent yield 
stress and lamina stiffness, respectively, and are to be determined 
experimentally. The yield stress, x^ , should be reasonably close to 
the matrix yield stress obtained from a test using matrix material 
alone as long as the damage occurs in the matrix rather than along 
the interface or within the fiber. The quantity Gj^/h is felt to be 
less well defined as discussed. 

The resulting nondimensional equations are: 

For all fibers, except N and N+1 when n , 


d^V. 


dn 


f ^ ^ = 0 


(B.19) 


for fiber N, n < « 


d^V. 


dn‘ 


- ''n + Vi -•^0'=''-^=' = ° 


(B.20) 
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and for fiber N+1, n ji a 



(B.21) 


d^V 


N+1 


dn 


■ ''n+1 ''n+2 + ° • 


Defining a new unknown function f(n) such that 

f(n) == V|^ ” - Tq < n-3 > if n < ot , 


(B.22) 


and 


f(ri) = 0 . 


ri > a 


with g(n) = '^N ■ ''n+ 1 for the same range of n values, the above three 
equations then become: 


d^V. 


+ v_., -2V„ + V. I = 0 , 


dn 


2 *n+l ‘■“n ’n-l 


d^V. 


dn 


2 * '^N+1 ■ ’ 


(B.23) 


(B.24) 


and 


d^V 


N+1 


dn 


2 ^^N+2 ■ ^'^N+1 '^N ^ 


(B.25) 


These differential -difference equations may be reduced to differen- 


tial equations by introducing the even valued transform as 

V 

0 


from which 


V(n»e) =^+ ^ V (n)cos(ne) , 
^ n=l ^ 


V„(n) = ~ / V(n,e)cos(ne)d 0 , 
n TT Q 


(B.26) 


(B.27) 


and the three equations become: 


f / f - 2[1 - cos( 9 )]v| cos(ne)d 9 = 0 
0 ^dn 


(B.28) 


for all fibers, except N and N + 1 when n < a , 
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(B.29) 


2 — ^ 

I* ~ cos(e)]V j cos(ne)cl 0 = - f(n) , 

0 V. dn > 

fiber N, n£ a, and 

2 — ^ 

— f^( - 2[1 - cos(e)]vl cos(ne)de = f(n) , (B.30) 

0 \ dr] J 

fiber N+1, n j< a . Making use of the orthogonality of the circular func- 
tions these three equations may be written as one equation, valid for 
all values of n and n , as follows: 

— f^( ^ - 2[1 - cos(e)]v) cos{ne)de 

"" 0 Idn^ J 

= — < a-n > f f (n){cos[(N + 1)0]- cos(Ne) } cos(ne)de . (B.31 ) 

^ 0 

This equation is of the form 

— /^F(Ti,6)cos(ne)de = 0 for all n and n (B.32) 

0 

and noting the definition of V(n»0) in equations (B.26) and (B.27) it 
is seen that the function F(n,0) is even valued in e and therefore, if 
the integral is to vanish for all n, the function F(n,0) must be zero. 

The single equation specifying V(nj0) is then 



- 6^V - - <a-n > D^f(n) 

(B.33) 

where 


6^ = 2[1 -cos(e)] = 4 sin'^(e/2) , 

(B.34) 

and 


= cos(Ne) - cos[(N + 1 )e] . 

(B.35) 


It is very significant that the irregular boundary condition, 
equation (B.l), of specified stress over a finite length, not coincident 


with either coordinate axis can be accounted for exactly and that the 
problem reduces to one differential equation which must satisfy boundary 
conditions along the coordinate axes only. The ability to do so strongly 
depends on the form of the failure criterion. A condition in which both 
normal and shear stresses were included generally would couple the axial 
and transverse equilibrium equations and yield a far more complicated 
set of differential equations. 

The solution to the problem of vanishing stresses and displace- 
ments at infinity and uniform compression on the ends of the broken 
fibers will now be sought. The complete solution is obtained by adding 
the results corresponding to uniform axial stress and no broken fibers 
to the following solution. 

The boundary conditions are then 

= 0 as n “ » (B.36) 

for all fibers. 


dV. 


= 0 . 


dn F|n 

for broken fibers, and 


= - 1 , 


for 


= 0 


(B.37) 


= 0 , for n = 0 , (B.38) 

for unbroken fibers. 

Using a technique such as variation of parameters to determine a 
particular solution to equation (B.33), the complete solution satisfying 
stresses and displacements at infinity is 


V(n,e) = A(e)e ^ < a-n >/°^sinh[6(n-t)]f{t)dt (B.39) 

n 

where the unknown functions are A(e) and f(t). The remaining two boundary 
conditions give 
20 
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cos(n0)de = - 1 


dV^(O) 



D^/ cosh(<5t)f(t)dt 


for all broken fibers, and 


(B.40) 


V (0) - — ( I\{q) - ^ /^sinh(6t)f (t)dtl cos(ne)d0 = 0 (B.41) 

for all unbroken fibers. Equation (B.41) is solved exactly by taking 

n2 a ^ 

A(e) - -T~ f sinh(6t)f{ t)dt = ^ B cos(me) (B.42) 

0 m=0 

where the B^ are constants. Equation (B.40) then gives a system of N + 1 
algebraic equations for the N+1 constants in terms of f(n) which 
is, as yet, unknown. For the case of no damage the problem is then 
exactly the same as section A of the present report. 

For matrix damage, a 0, equation (B.40) must be supplemented by 
the condition that 

f(n) = g(n) " ^0 ^ ^ ’ n £ a 


= < n-3 > , (B.43) 

and recall from equation (B.22) that f(a) = 0 and therefore 

g(a) = Tq . (B.44) 

The constants B^^ and the function g(ri) are then specified by 
requiring that equations (B.40), (B.43) and (B.44) be satisfied. Using 
equation (B.39) and after considerable algebraic manipulation, the 
displacement of any fiber for all values of n is 
P ^ . N 

V = — / e z cos(m 0 )cos(ne)de 
" 0 m=0 

+ 7 A(t){C^(|t-nl) -C„(t+n)}dt , 


(B.45) 



where 


C^(C) = f ^ cos(ne)de 
Equation (B.40) then becomes 

~/f-6 z B cos(m0) + D^/“e"'^^g(t)dt-D^ t /“e'^^^dt 
^ 0 I m=0 0 ° 3 

X cos(ne)d0 = - 1, n = 0,l,...,N (B.46) 

and equation (B.43) along with (B.45) gives, for n £ ct , 
g(n) = , 

= — z B cos(me){cos(N0) - cos[(N + 1 )0]}d0 

"^0 m=0 

+ Y f g(i') “ C^(t+Ti) “ ( I ! )"^ ^N+1 

" ^ ^ ^ ^ “ ^N+i ^ I I ^ ^ ^N+1 * 

(B.47) 

The condition that 

g(a) = Tq (B.48) 

also must be satisfied. 

Physically, it would be more direct to specify the applied stress 
and the number of broken fibers, N, and determine the damage zone 
a and 3 depending on given yielding and splitting conditions. As a and 3 
appear in the limits of the above integrals this is not convenient 
mathematically and it is easier to specify the number of broken fibers, N, 
and the damage zone a and 3» and compute the required applied stress . 
These equations were solved as follows: 
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(I) An initial set of constants was determined for the problem 

of no damage, a - 3 = 0 in equation (B.46), i.e., 

N p TT 

z B — / 6 cos(me)cos(ne)d 0 = 1 , n = 0,l,...,N. (B.49) 

m=0 0 

(II) These initial constants were then substituted into the integral 
equation (B.47) and, along with equation (B.48), the function 
g(n) and x^were determined using the desired values for a and 3 . 

(Ill) Using g(n) and x^., a new set of constants, B^ , was computed from 
equation (B,46) with the desired values of a and 3 • 

(IV) This procedure was repeated until the unknowns changed less 
than a prescribed amount with additional iterations. 

In the above solution the unknown function, g(ri), was assumed to 
be piecewise linear over the interval 0 £ n £ a of the form 

(n) = + yj n , i = 1 ,2, . . . ,k 

when the interval was divided into k equal subdivisions. The function 
g(n) then contained 2k unknowns with one additional unknown being Xq . 

As g(n) is the displacement difference it should be a positive, mono- 
tonically decreasing function and its representation as a piecewise 
linear function should be sufficiently accurate. The ( k+1) equations 
were obtained by requiring that the integral equation, equation (B.47), 
be satisfied at the (k + 1) end points, (k-1) equations resulted from the 
requirement of continuity of the function g(n) between adjacent interval 
and the last equation was given by g(a) = x^. 

With the longitudinal displacement v^ now known, the transverse 
displacement u^ is obtained by solving equations (B.4), (B.6), and (B.8). 
This solution is recorded below for completeness. 



sin(n9 )d9 
1 - cos(9) * 

(B.50) 

The matrix normal stress, aw! > is related to the transverse displace- 

M| n 

ment, , through the stress-displacement relations discussed in sec- 
tion A. Knowing u^ from equation (B.50) the matrix normal stress can 
be computed for all values of n and n. 

In the next section this solution is modified to account for 
transverse damage. 


“n = - 


o G|^h 

'n’EniiEr 
M F 


- 

^ ■ 2"^ {sin(Ne) + sin[(N + l) 0 ]} 
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^ • Two-Dimensional Shear- Lag Model with Broken Fibers, Longitudinal 
r i x_ Spl itting and Vieldinq and Transverse Matrix and Fiber 
Damage 

A certain amount of stable transverse extension of the initial 
notch under increasing applied loads has been observed in tests on 
unidirectional and cross-ply composite laminates [13]. In case of 
single-ply and multi-ply unidirectional laminate this occurs in the 
form of breakage of an arbitrary number of fibers ahead of the initial 
notch tip sometimes accompanied by fracturing of the matrix and more 
often without any matrix fracture. In the later case the damage 
may be observed with x-ray or by etching away the matrix. The matrix 
material undergoes extensive longitudinal yielding in both the cases. 
This behavior appears to be strongly dependent on the laminate thick- 
ness; a detailed investigation into this question is currently underway 
by the present writers. The breakage of fibers is found to not be 
confined to a straight line in the transverse direction, that is, all 
fibers do not break along the x-axis even though the initial notch 
is oriented along the x-axis. This creates a zone ahead of the 
notch tip in which fibers do not possess their original stiffness and 

hence results in a reduced axial load carrying capacity of these 

fibers. If all the fiber breaks are not along the same line this 

reduction will not be as drastic as for rectilinear extension and 

thus the damage zone will support considerably more load than that 
of a matrix region alone. 

In order to account for this behavior, an idealized model having 
transverse damage ahead of the initial notch in addition to the longi- 
tudinal matrix damage considered in section B is now developed. 



This is shown in Figure (4). It is mathematically untractable to 
account for the distribution of breakage of fibers occurring at 
points other than those on the x-axis; the model then assumes all 
breaks to occur on the x-axis and accounts for the stiffness by 
assuming these fibers to carry a reduced load as described below. 

The transverse damage consists of an arbitrary number of broken fibers 
which are constrained by the adjoining matrix and/or by the unbroken 
fibers through the thickness. These fibers in the transverse damage 
zone will be referred to as constrained fibers. The extent of con- 
straint is represented by a stiffness coefficient y and is assumed 
to be constant for all the constrained fibers. The stiffness coeffi- 
cient, Y , is given by 

stress in the constrained fiber 
^ ~ stress in the first unbroken fiber 

Longitudinal damage, yielding and splitting, is assumed to occur at 
the end of the original notch as in section A. 

With reference to Figure (4), n = 0,l,...,N corresponds to broken 
fibers in the initial notch region, n = N + l, N + 2,...,M corresponds to 
constrained fibers in the transverse damage zone and n = M+l, M+2,...,°° 
corresponds to unbroken fibers. By comparing Figures (2) and (4), it 
may be observed that the only difference between this model and the 
one considered in section B is in the boundary conditions along the 
x-axis. These boundary conditions are given in Figure (5). The 
governing differential -difference equations and also the final equa- 
tion specifying will be the same as those of section B. Starting 

from the differential equation derived in section B, the solution is 
obtained using the appropriate boundary conditions following the pro- 
cedure described in section B. 
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The governing differential equation is given by 


071 “^ 

where all the quantities have the same meaning as in section B. VJith 
reference to Figure (5), the boundary conditions, after normalization 
are given as follows: 


(i) 

As n ^ 

(ii) 

At n “ 0 


and 


= 0 for all fibers. 


dV 

= - 1 for broken fibers, 

n 



(C.2) 

(C.3) 


dV 
n 

dn 


-an - - 1 + Y On 
'■ n r 


M+1 


for broken fibers, (C. 4) 


where a, 


is the normalized stress in the first unbroken fiber at 


f|m+i 

ri = 0 and Y is the stiffness coefficient defined earlier. 


(iii) At n = 0 = 0 


for unbroken fibers. 


(C.5) 


As in section B, the complete solution to equation (C.l) satisfying 
vanishing stresses and displacements at infinity is given by 

2 

V(n»0) = A(e)e~'^^ + ^ <a-3 > sinh[6(n-t)]f (t)dt , (C.6) 

n 

where the unknown functions are A(e) and f(t). The remaining two bound- 
ary conditions give 


dVn(O) ^ 2 i'. „ 

' IT Q |-6A(e) + D^ cosh(6t)f(t)dtlcos(n0)d0 = -l 

n = N+l,...,N, (C.7) 


dn 


dV^(O) 

dn 


2 ^TT 
TT 


f /-6A(0) +D^/ cosh( 6t)f ( t)dt \ cos(n0)d0 = -l +y5_ 

0 I 0 J 


M+1 

n=N + l,...,M, (C.8) 
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and 


-.(0) . I /{, 


A(0) - 


^ / sinh(st)f(t) } cos(ne)d 0 = 0 


n = M + l, M + 2,...,a> 


Equation (C.9) is solved exactly by taking 


(c.9) 


2 M 

A(e) - ^ /°*sinh(6t)f(t)dt = z cos(me) 

^ 0 tn=0 

where are constants. The stress in the first intact fiber. 


M+1 


, is given by 


(C.10) 




1+1 9ot I 

^ ^ cosh(6t)f(t)dt j cos[(M+l)e] 


= - 1 + 


M+l 


(C.llj 


Using equations (C.IO) and (C.ll) in equations (C.7) and (C.8), A(e) 


and Oi 


M+1 


may be eliminated resulting in M + 1 algebraic equations for 


M + 1 constants in terms of f(ri) which is, as yet, unknown. For 
longitudinal matrix damage, equations (C.7) and (C.8) must be supple- 
mented by the conditions that 

f(n) = g(n) - < n-3 > > n < a , (C.12) 

and 


g(a) = 


(C.13) 


From equations (C.6) and (C.10), A(e) may be eliminated to obtain 
V(rij0) in terms of constants and unknown function f(t). Recalling 
the relation between V(n>0) and V^(n), an expression can be obtained 
for the axial fiber displacement V^(n) as 

p TT _(5r| 

V„(n) = — / e 2 : B cos(me)cos(n0)de 

^ ^0 m=0 ^ 


+ i/f(t){C„Ut-Til) -c lt+Ti)}dt , (C.14) 


28 



where 


C„(e) = - cos(n6jd0 

n 0 ° 


Equations (C.7) and (C.8) then become 


(C.15) 


s B cos(me)+D^/ e g(t)dt-D^x / e”'^^dtl 
0 V m=0 0 ^3 J 


cos(ne)d 0 = - 1 


n = 0,1,. ..,N, (C.16) 


and 


9 TT 

■ i 


M 9 ot -cSt _ 9 -6t 

-6 s cos{me) + D / e g(t)dt - t D / e dt 

m=0 "" 0 ° 3 

{cos(ne) - Y cos[(M + l)0]}de=-l+Y n = N + l,...,M. (C.17 ) 


Equation (C.12) along with equation (C.14) gives 

M 

q(n) = — 2 B„ cos(m0){cos(ne) - cos[(N+l )0]}d0 

0 m=0 

+ Y |C|^(|t-nl) -C,^(t+Ti) - |t-n|) + (t+n)j dt 

- ^ /“fc^{|t-nl) -Cf,(t+n) -C|^+l(|t-n|) + Cfj^.i(t+n)jdt. (C.18) 
The condition that 

g(a) = (C.19) 

must also be satisfied. 

Equations (C.16) to (C.19) are of the same form as those obtained 
in section B for g(n) and the constants B^ ; however, they now contain 
the additional parameters M and y . The procedure to obtain displace- 
ments and stresses is the same as before. For boron/aluminum laminate 
the size of the damage zone (M-N) was found to be approximately constant 
for all initial notch lengths (N). The coefficient y decreased with an 



increase in the initial notch length. The results are d 
where it is shown that table transverse notch extension 
with the extension becoming unstable in a boron/aluminum 
about seven damaged fibers. 


scussed below 
s possible, 
laminate at 
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D . Two-Dimensional Shear-Lag Mod e l with B ro ken F i bers. Longitudinal 

Matrix Splitting and Yielding with a Surface "Cbnstraint Lay^ 
(Cover Sheet) 


A constraint layer is added parallel to the main lamina of the 
problem considered in Section B, providing additional shear stiffness. 
The constraint layer is an attempt to account for either the misalign- 
ment of fibers in a multi -ply unidirectional laminate, or the presence 
of a non-zero ply in the composite laminate. The extent of the above 
effects may be represented by the amount of constraint the layer intro- 
duces, defined in Figure (6) as the constraint ratio. The constraint 
ratio may be varied by varying the constraint layer parameters. How- 
ever, in the present investigation no explicit relationship between 
the constraint ratio and constraint layer parameters is considered. 
Displacement-stress relations similar to those in section B are assumed 
so as to obtain decoupling of the equilibrium equations. Only fiber 
stresses and matrix shear stresses will be developed in the present 
study. The equilibrium equations in the transverse direction are more 
difficult than those in the previous models due to the presence of the 
constraint layer and will be considered at a later time. 

A free-body diagram of an element consisting of the main lamina 
and the constraint layer is shown in Figure (6). With reference to the 
free-body diagram, the equilibrium equations in the axial directions are 
given below. For all fibers n, with the exception of N and N+1 when 
y ;< L the equilibrium equation is 


Ac dCTn „ 

F F n , 

T“a3r~ ■^'^i 


- T 


= 0 . 


n+1 


(D.l) 


For fiber N and N+1 , y L, t 
tions are 


N+1 


= - Tq <y-z> and the equilibrium equa- 
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t dy 


•|n 

- T < y -^ > - T 


1 “ T ' 

y O 

N 1 

In+1 




(D.2) 


and 


Ap dap 

— _^+ + Tq <y-)l> + -T 


N+2 


t 


N+2 


N+1 


}^= 0 


(D.3) 


The following stress-displacement relations are introduced in equations 
(D.l) through (D.3) as before: 
dv. 


= E. 


n 


Tin F dy ’ 

Gw 

^In+l =TT(Vi - ''n)’ a"'' 


(D.4) 


ln+1 h' ^''n+l " '^n^ ‘ 

The equations (D.l) through (D.3) then become 

ApEp d V Gw p'f' 

+ Tr(Vl-2''n + Vl) W(Vl-2Vn + v„.i) = 0 (D.5) 


' dy2 


for all fibers, except N and N+1 when y ±l , 


ApEp d V 1^ Gw p I ^ I 

■ ■'o ■ T ^''n ‘ ''N-1^ ^''n+1 ■ ^''n ''n-1 ^ ° 

dy 


for fiber N when y £ L, and 


(D.6) 


ApEp d 


6, 


dy 


2 ' + ^0 < y-ii > + h " '^N+l ^ W ^''n+2 ■ ^''n+I ''n^ " ^ 


(D.7) 


for fiber N+1 when y ^ L . 


G ^ t ^ h 

By defining ^ T~ ^ r constraint ratio, and rearranging terms 

equations (D.5) through (D.7) may be rewritten as 
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and 


ApEph 

d^v 

+ 

V 

dy^ 

ApEph 



? + 
dy^ 

ApEph 

‘^^''n+1 


dy^ 


(1 + - 2v„ + v„.i) = 0 , (D.8) 

(1 + CR)(Vf,^.i - 2 v,^ + Vf,_i) = ^^0 <y-^> + Vf(+ 1 -Vfj . 

(D.9) 

+ (1 +Cp)(Vf,^2"2vN+i +V n) = ''n"''n+1‘ ^'^0 • 


(D.IO) 


The same change of variables as in Section B will be used in order to 
normalize the quantities in the equilibrium equations. The resulting 
nondimensional equilibrium equations are: 


for all fibers, except N and N+1 when n a , 
d^V 

-2V„ + V„.,) = 0, (D.ll) 

dn 

for fiber N, when n ^ a 


2^ + (1 + Cf^)(V^^-j - To<ti-3>, 

dn 

and for fiber N+1, n ji a 

2 n + CR)(Vf|^,2 - ■ ■ ^'^N+1 ■ ‘ 

dn 


(D.12) 


(D.13) 


Defining a new function f(n) such that 


f (n) 

and f(n) 

with g(n) 
equations 


= Vj^ - V|^^i - TQ<n-3> if n < ct 


= 0 


if n > ot , 


= Vn - , for the same range of n values, 

then become: 


(D.14) 
(D.15) 
the above three 
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(D.16) 


d^V 


dn‘ 


r + (1 + Cn)(V 


R'''n+1 


d^V, 


dn 


f + (1 + Cj^)(V 


N+1 


and 


d^V 


dn 


^ + (1 + Cp)(Vm^2 






2Vn+1 + V,) 


= 0 , 


= - f(n). 


(D.17) 


= f(n) 


(D.18) 


These differential -difference equations may be reduced to differen- 
tial equations by introducing a new function V(n»0) previously defined 
in section B, so that, Vn(n) is given by 


Vp(n) = - / V(n,e) cos(ne)d0 . 
^ 0 


(D.19) 


After substituting for V^(ri) from equation (D.19) the above three equa- 
tions become, respectively. 


/ < — 5 - 2(1 + Cd)[ 1 - cos(e)]V > cos(ne)de = 0 , 
0 I dn^ ) 


(D.20) 


2 j.'' r d^V 

7T 


and 


0 Kdr] 


^ d^V 


2 - 2(1 + C|^)[l - cos(e)]V } cos(N 0 )de = - f(n) , 


"" 0 ^d 


/ < ^ - 2(1 + C|:^)[l - cos(0)]V > cos[(N+l )n]d = f(n) 


(D.21) 


(D.22) 


Making use of the orthogonal i ty of the circular functions these 
three equations may be written as one equation, valid for all values of 
n and n as follows: 

- f - 2(1 + Cn)[l - cos(o)]V^ cos(n0)d0 

0 Idn ^ J 

= — < a-n > / f (n){cos[(N+l )0] - cos(N0) } cos(n0)d0 . (D.23) 
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This equation is of the form 


2 

— r 

^ 0 ^ Q » '^0'^ 3ll r) and n 

Since F(n,0) is even valued in 0, if the integral is to vanish for all 
n the function, F(n>0) > must be zero. The single equation specifying 
V(ri 9 0) is then given by 


_ 6^ V = - < a-n > D^f(n) , 

(D.24) 

dn^ 

= 2(1 +C|^)[1 - cos(e)]= [2 / 1 +Cp s1n(e/2)]^ , 

(D.25) 

= cos(Ne) - cos[(N+l )e] 

(D.26) 


The equation (D.24) is exactly of the same form as the correspond- 
ing differential equation, equation (B.33) of section B, the only dif- 
ference being in the definition of 6. Consequently, the expression 
V^(n)» the algebraic equations for the constant B^, and the integral 
equation are the same as in section B. They are reproduced here for 
completeness. 

The displacement of any fiber for all values of n is 

Vj^(n) = — / e~ z B|^ cos(m0)cos(n0)d0 
0 m=0 

+ J f(t)|c^{ |t-n|) -C^(t+Ti)jdt , (D.27) 

where 

c„(c) = f Z X e’'^'cos(ne)de . 

The algebraic equations for the constants B^^ are given by 



2 TT^ M 
- / <- 6 E 
0 I m={ 


. 9 cc _j-|- o ^ oc _P4- 

B cos(me) + D / e g(t)dt - D t / e ^ dt? x 
0 

cos(n0)dt = - 1 , n = 0,1,.. .,N . (D.28) 


The function g(n) is given by the integral equation 


g(n) = — / E cos(m9)D^de 


+ ^ / g(t) |C(^( It-nl )-C|^(t+n)-Cj^^^ ( I t-n| )+Cj^^-j (t+n)| dt 

- “2“ ^ “ ^|^(t+n) “ ^^+'1 ( 1 1 dt . (D. 29) 

The condition that 

g(a) = » 

must also be satisfied. 

The numerical technique described in section B is employed to 
obtain the solution to the above equations. 
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E. Two-Dimensio na l Shear-Lag Model with Broken Fibers, Longitudinal 
Matrix Split ti ng, and Yielding w ith an Elastic-Linear Strain- 
Hardening Matrix 


In section B, longitudinal matrix splitting and yielding was con- 
sidered with the matrix material being elastic-perfectly plastic. In 
this section a bilinear stress-strain behavior for the matrix material 
is assumed, as shown in the figure below. 



Bilinear stress-strain relationship 


A free-body diagram for a typical element is given in Figure (2) 
with the special condition for the last broken fiber, N now given by 


T 


N+1 


= CG*Yf^+1 


(^M - G*)yq] <y-ii > 


(E.l) 


where the negative sign is taken to account for negative shear strain, 
so that the absolute value of Yq is specified. 

The egui librium equations in the longitudinal and transverse 
directions, respectively, for all fibers n, with the exception of N 



Ar dap „ 

■f ■^In+l ^In = ° ’ 


h d 


'M n+1 


■ '^H n 2 dy ^'^ln+1 ^ ’^In' 


= 0 . 


(E.2) 


For fiber N, y < L the equilibrium equations are 


Ap dap I 


f- [G*y„,i - (G„-G*)y,] <y-.> ' t|^ = 0 . 


+ ]=0. (E.3) 


For fiber N+1 , y £L the equilibrium equations are 


Ar- da 




°m1n+2’^m|n+ 1 2 dyl^^L+2 

(E.4) 

Introducing the stress-displacement relations defined in sections 
A and B, and also noting that = (v^^.-j -v^)/h, the equilibrium 
equations in the longitudinal direction may be written as follows. 

For all fibers except N and N+1 when y :i L , 


ApEph dS„ 

^ ( Vl - + v„.p =0 . 

For fiber N, when y £ L , 


G^t N I N bpi 0 


Gm Gpi 0 


G* 


(E.5) 


(E.6) 
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For fiber N+1 , when y l. L 


ApEph d 

,2 ''n+2 ■ ''n+1 




dy‘ 


+ (1 - ^ T <y-£> 

'"M “ 

■ It ^''n+1 ■ ''n^ < y-s. > = 0 


(E.7) 


The equilibrium equations in the longitudinal direction are independent 
of transverse displacements u^. Therefore, only the solution for the 
longitudinal displacement will be considered. 

The same change of variables as in section B will be used in 
order to normalize the quantities in the equilibrium equations. The 
resulting nondimensional equations are: 
for all fibers, except N and N+1 when n j< a 


d^V. 


dn 


f " ''n+l - = 0 


(E.8) 


for fiber N, when n < « 


d^V 

'^N-l " ('^N+r '^N^ < n-3 > = 0 , (E.9) 


for fiber N+1, when n < a 


2 

V 

~‘2^ '^N+2 ■ Vl ■^o ^ ^ ^''n+ 1 - Vfj) < n-B > = 0 . 

(E.IO) 

p'k 

For the sake of simplicity, let Gq = ^ and 

M 


- n - 
'"or '"o • 

As before, a new function f(n) is defined such that 

f(n) == Vj^ - - Gfj(Vj^+-j - V^)j < n 6 > , if n < a , 

and f(n) = 0 , if n ^ a . 


(E.ll) 


(E.12) 
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With the introduction of f(n), the equations (E.8) to (E.IO) may 
be written as follows: 




dn 


Vl-^Vn + Vl =0 


d^V, 


and 


dn 


d^V 


2 '^N+1 ■ ’ 


N+1 


dti 


2 ''n+2 ■ '^N • 


(E.13) 

(E.14) 

(E.15) 


Equations (E.13) to (E.15) are exactly of the same form as the cor- 
responding equations of section B, i.e., equations (B.23) to (B.25). 

The only difference is in the function f(n)» which is now defined as 

f(n) = g(n) - < n-3 > » for n < a , (E.16) 

so that 

g(n) = (1 -6j^ Vl^ > 


and 


g(a) = 


(E.18) 


oR * 

The boundary conditions on stress and displacements are the same 
as those of section B. Following the procedure detailed in section B, 
the axial displacement in fiber n for all values of n is given by 

O 7T — n ^ 

V (n) = — / e ^ z B cos(me)cos(n9)de 
^ ^ 0 m=0 ^ 


J- /*f(t)|c^( |t-n|) - c^(t+n)jdt , 


(E.19) 


where 




Cn(^) = f Y ^ ’ cos(ne)d0 , 

6=2 sin(0/2), 

and q 2 „ cos(Ne) - cos[(N+l)e] . 
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Boundary conditions on fiber stress yield the following equation: 


o 7T r ^ ? a -54- p a -.fit 

-fj-6Z cos(me) + D / e ^g(t)dt-D f^n / e dt)X 

^ 0 V m=0 0 B J 

cos(ne)de = - 1 , n = 0,1, ...,N . (E.20) 

Putting n = N and N+1 in (E.19) and subtracting from V|^ gives 

p TT N p 

+ Y /“ |t-n|) -c^(t+n) -C„+i( It-nl) +c^^.^(t+n)|dt. 

(E.21) 

Multiplying equation (E.21) by (1 - < n -3 > ) and recalling the defini- 

tion of g(n)» equation (E.21) reduces to the following: 

p T7 -fin N p 

g(n) = [1 - Gn < n-3 >] - / e ^ B„ cos(me)D de 
^ l'^ 0 m=0 

+ j n-3 >j It-n 1 ) - Cj^(t+n) - ( It-nl) 

+ C|^(t+n)^dt . (E.22) 

Further simplification of (E.22) yields the integral equation 
g(n) = [l-GD<n-3>] - sB^ cos(me)cos(ne)D^de 

K "^0 ^ 

+ ^ / g(n)^C^( It-n l) - Cf^(t+n) - Gj^j+i ( |t-nl) + ( t+n)| dt 

- -p- f ^Gj^( 1 t-n 1 ) - G^( t+n) - ( It-nl) + ( t+n)^ dt 

(E.23) 

The last condition that must be satisfied is 
9^“^ ” ^OR ^ "^0^^ " * 
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Equations (E.20) and (E.23) together with equation (E.24) are solved 
employing the numerical technique described in section B in order to 
obtain displacements and stresses. 

The above development of the different modifications to the fun- 
damental solution of section B has been presented with each change con- 
sidered separately in the past three sections (C,D and E). This pro- 
cedure was used to note more clearly the necessary changes. However, 
in developing the computer codes for these solutions we have included 
all the modifications in one code and have investigated the influence 
of the various models as a complete set. In the results, to be dis- 
cussed later in this report, it was found that the transverse damage 
model of section C was the most significant. The inclusion of a strain- 
hardening matrix had very little influence on the ability of the solu- 
tion to give results consistent with experiments for ductile matrix 
laminates. Perhaps if unloading had been considered, more dependence 
would have been noted. Also, the constraint layer (cover sheet) was 
not significant in comparing with experimental data on unidirectional 
laminates, however it certainly would be for cross-ply laminates. 

The following two sections present the solutions for a rectangular and 
circular damage region without including a strain-hardening matrix or 
a cover layer but including all other damage. The main purpose was to 
investigate the changes given by different initial damage shapes. 
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F. Two-Dimensional Sh ear-Lag Model with a Rectangular Notch, Longi- 
tlJdina l' M atrix SplU^ and Yielding and Transverse Matrix and 
Fiber Damage 

In the next two sections we extend the above solutions to account 
for an initial damage zone in the shape of respectively ; a rectangular 
notch and a circular cut-out. These represent more realistic damage 
shapes and it is of primary interest to investigate the differences 
in the results as compared to the idealized model of a slit as previously 
discussed. The same basic assumptions as before will be made and, as 
the solutions are developed very easily from the past work, only the 
fundamental differences will be discussed. Both longitudinal and trans- 
verse damage will be included initially. That is, the development will 
concentrate on the necessary changes to the complete solution presented 
in section C rather than starting with no damage and developing each 
successive solution. For both of these cases the only significant 
difference in the solutions is the change in boundary conditions on the 
initial damage region and, for the rectangular opening in particular, 
these differences give almost trivial changes in the resulting form 
of the equations. 

The solution to be developed in this section is for an initial 
damage region having the form of a rectangular opening. Under loading, 
damage is taken, as before, in the form of longitudinal splitting and 
yielding of the matrix and transverse damage to the fibers. The longi- 
tudinal damage is assumed to occur at the end of the initial notch 
between the last broken fiber and the first unbroken fiber. Transverse 
damage is modeled, in the same manner as in the preceding sections, as 



reduced load carrying fibers along the horizontal axis. Figure (7) 
depicts the geometry for this study. This last assumption is admittedly 
a simplification as one would assume that the location of the maximum 
axial stress in the first unbroken fiber in front of the notch, before 
transverse damage occurs, would be at the corner of the notch rather 
than the center. If this is true then notch extension would, at least, 
begin at the corner. The results from the present solution indeed in- 
dicate the maximum stress to occur at the corner. However, the difference 
between this maximum stress and the stress at y=0 is small and further, 
the location of the maximum stress in adjacent unbroken fibers away 
from the notch changes from the corner to the y-axis after only a few 
fibers. Preliminary experimental studies now being initiated also 
indicate that transverse extension does originate at the notch corners 
but successive breaks occur in a random manner tending to be symmetric 
about the x-axis. So, even though the model does not account for this 
irregularity in transverse extension, it is felt to be a reasonable 
approximation, especially for larger notch extension. 

With the above assumptions, the solution is then developed by 
modifying the previous results to account for the boundary conditions 
of zero shear stress on the sides of the notch and a unit compression 
stress on the ends of the fiber at the top (and bottom) of the notch. 

As the solution already allows for matrix splitting, the first of the 
above conditions is met simply by setting the minimum value of the split 
length ( 3 ) equal to the half-notch height. The remaining boundary 
condition is satisfied by taking n = H (normalized half-notch height) 
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^ |ro 


rather than n=0 as in equation (B.46J; this is 


IT N _ pU 

/ - E B cos(m0)5e 
0 m=0 


cos(n 0 )de = - 1 


All other equations are unchanged. 


G . Two-Dimensional Shear-Lag Model with a Circular Cut- ou t, Longi- 
tudinal Matrix Splitting and Yielding and Transverse Matrix 
and Fiber Damage 

The solution to this problem for the case of no additional 
damage was presented by Franklin in [7] and Kulkarni et. al . in [8]. 

Both of these studies had difficulty satisfying the stress free 
boundary condition of the circular boundary; Franklin's solution gave 
zero stress in the fibers but did not consider the non-zero shear 
stresses in the matrix while Kulkarni attempted to remove the shear 
stresses by an averaging technique. The results of the two are, 
however, not significantly different. In the present solution we 
formulate the boundary conditions in the same way as Franklin but 
the ability of the present model to allow for matrix splitting gives 
some necessary additional freedom to satisfy the stress free conditions 
more accurately. Longitudinal and transverse damage is also included 
with the longitudinal damage being at the edge of the hole and the 
transverse damage originating and extending along the horizontal 
axis as in the slit and the rectangular notch. 

To develop the stress free boundary conditions on the hole 
surface first consider an unnotched laminate having uniform applied 
stresses at infinity and determine the fiber stresses on a circular 
region having the radius of the desired cut-out as shown in Figure (8). 
The negative of these stresses are then the appropriate boundary 
stresses to be applied to the edge of the circular hole in an infinite 
laminate having vanishing stresses and displacements at infinity. 

A typical element on the boundary of the circular hole is shown 
in Figure (9). A series of these inclined elements were joined together 
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to form an approximation of the circular boundary as shown in Figure 
(10). It is seen in Figure (10) that for each element a small vertical 
matrix boundary exists on which the shear stress, as given by the shear- 
lag assumption, will not be zero. The largest such surface is between 
the last broken and first unbroken fiber. By setting the matrix 
split length ( 3 ) equal to this distance, this portion of the shear 
stress boundary condition is satisfied. The shear stress condition 
on the remaining elements is not satisfied. 

Referring to Figure (9) and summing forces in the horizontal and 
vertical directions yields 

-^Vn -t^RT = 0 , (G.l) 

and 

-tVn cos(8n) + tT^-r sin(0j = 0 . (G.2) 

Solving for the boundary stresses yields 

Af 2 

°R = Ht si" (%) » (G.3) 

and 

ht si"(en)cos(0j- (G.4) 

As previously mentioned, the negative of these stresses will be 
applied to the boundary of the hole in order to remove the resultant 
stress on the hole boundary. Since it is only necessary to solve the 
axial equation in order to obtain the axial stresses and displacements, 
a boundary condition in the axial (fiber) direction is the only one 
needed. Thus the appropriate boundary condition can be derived by 
summing forces in the axial (y) direction. 


ht 


sin(e^) + T 


n+1 


cos(e^) = 


- Rt 


(G.5) 


n 


47 



The same equation can be derived by summing forces in the radial and 
transverse directions and solving the resulting equation simultaneously. 

The stresses can now be expressed in terms of displacements by 
using the shear-lag assumptions as before. Upon substitution, the new 
boundary condition becomes 

ApEp dV -a Ap 

T5j;jr -3F tVl -V sin(e^) . (G.6) 

From the geometry of the hole, as seen in figure (8), the radius 
can be expressed as 

R = h, (G.7) 

where the coordinates of a particular point are 
x^ = hn and 

= h _ (G.8) 


The boundary condition for the circular hole then becomes 


dn 


a. 




where 


2n 

% 2N+1 


and 


[V 


n+1 


-V = - 1 





(G.9) 


This boundary condition must be satisfied at the appropriate value 
along the boundary of the hole. It is seen that the normalized coordinate 
n is then 


n(nj 


'^n _ 2N+1 r; 2 

hH^ 2H^ \ ' ■ “n 
0 0 * 


(G.IO) 


48 



Finally, the boundary condition becomes 


2 IT N 

— / - D B cos(m0)(5e cos(n0)de 

7T ^ r\ III 

0 m=0 

2 TT ^ 

— / E B cos(m0}e ^{cos[n+l )0] - cos(n 0 ) } = - 1 . 

"" 0 m=0 ^ 

(G.ll) 

This boundary condition replaces equation (B.46) used in the slit prob- 
lem. Note that this solution is not independent of material proper- 
ties due to the term defined above. 

Since the damage is defined in the same manner as in section B for 
the slit, the integral equation for the present solution remains the 
same. The only other change is that the slit length (3) must be, at 
least, equal to the matrix mismatch on the last inclined element 
This value is given by 


- 

“1 


^-“n_ 



The solution is obtained in the same manner as in the previous sections. 



H . Two-Dimensional Shear-Lag Model with Broken Fibers for the Half-Plan e 


The motivation for much of the above work has been the interest 
in attempting to develop approximate analyses capable of predicting 
accurately the fracture behavior of hybrid {buffer strip) composite 
laminates. In this report solutions have been developed and the sig- 
nificance and validity of the various models will be demonstrated for 
the unidirectional full -plane laminate following this section. Based 
on the understanding developed in this phase of the study, future work 
will involve the extension of these models to mul ti -material laminates. 

As a first step in this direction the solution for adjoining uni- 
directional half-planes containing damage at, or near, the interface 
will be investigated. The solution presented below is a fundamental 
part of this extension, and consists of the solution for an edge crack 
in a unidirectional half-plane. Using superposition, i.e.; matching 
boundary conditions on the interface, this solution can then be added 
to a second half-plane having different fiber and matrix properties to 
construct the adjacent half-plane problem. 

A two-dimensional array of parallel fibers with an arbitrary num- 
ber of broken fibers at the free edge is shown in Figure (11). The 
laminate is subjected to a prescribed shear stress and transverse normal 
stress along the free edge in addition to a remote uniform tensile stress 
in the axial direction. No additional damage other than the broken 
fibers at the free edge is considered at this time. 

With reference to the free-body diagrams shown in Figure (12) the 
equilibrium equations in the axial and transverse directions are given 
below. The boundary conditions on the free edge must now be specified 
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as opposed to the full -plane solution, in which they were accounted for 
by symmetry. We assume the first fiber in the half-plane to be embedded 
in the matrix material and, as seen in Figure (11), it will have a 
different set of equilibrium equations than all the remaining fibers. 

The equations of equilibrium are: 



(H.l) 

(H.2) 


for the first fiber. 


and 


Apdop 


+ t T 


n+1 


- t 


= 0 


n+1 


- a, 




(H.3) 


(H.4) 


for all other fibers. 

Substituting the stress-displacement relations of Section A into 
the above equilibrium equations the following pairs of equations are 
obtained : 


ApEph dS. . 

(V2-V^)-Tg ^ 


Gi^ t 

M 


^ = 0 , 


dy‘ 


M 


and 


E G 

ir(u9“U-j)-ap^|l+-^ 7 ^(v9-Vt) + y7T;7'r;, = Os 


dy 


1 


2 dy 


(H.5) 


(H.6) 


for the first fiber. 


ApEph dS^ 




dy‘ 


(H.7) 


and 
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(u 


n+1 


“ 2u + u 1 ) 
n n-r 




n+1 


- ''n-1 ) = 


(H.8) 


for all other fibers. 

From equations (H.5) and (H.7) it can be observed that the axial 
equilibrium equations are independent of transverse displacement. 
Since only fiber axial displacements and axial stresses will be 
developed in this section, the solution for the transverse direction 
equilibrium equation will not be given here. 

The same change of variables as in section A will be used to 
normalize the quantities in the equilibrium equations. The resulting 
equations in the nondimens ional form are given by 


d^V. 


dn' 


+ (V,-V,) - ^ = 0 , 


for the first fiber, and 
dn 

for all other fibers, where 


(H.9) 






Defining a new function f-|(n) such that 

n(n) = . 

for all values of n the above two equations become 


d^V. 


V 2 - 2 V, = f^(n) 


for the first fiber, and 


(H.IO) 


(H.n) 
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(H.12) 


d^V. 


dn 


Vl-2Vn + Vl = 0 


for all other fibers. 


These differential -difference equations may be reduced to dif- 
ferential equations by introducing a new function V(n»0) such that 
the normalized displacement V^(ri) is the Fourier coefficient in a 
Fourier series expansion. In order to represent the free-edge boundary 
conditions as discussed earlier, it is appropriate to utilize a sine 
transform. Then 


and 


V(n»e) = ^ Vj^(n)sin(ne) , 

n=l 


(H.13) 


V„(n) = - V(n,9)sin(ne)de 

7T Q 


(H.14) 


With the introduction of V(n>e) the equilibrium equations become, 
respectively, 


- 2[1 - cos(e )]v| sin(e)de = f-|(n) 


(H.15) 


and 


|- / / - 2[1 - cos(e)]V I sin(ne)d0 = 0 

0 t dn 


(H.16) 


Making use of the orthogonality of circular functions these two 
equations may be written as one equation valid for all n as follows: 


2 d^V 


I 


dn‘ 


2[1 - cos(e )]V 


sin(ne )de 

= — f,(ri)sin(e)sin(n0)d6 . 

TT Q I 


This equation is of the form 
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for all Ti and n , 


— F(n,e)sin(n0)de = 0 , 

0 

and noting the definition of V(n,0) in equations (H.13) and (H.14) it 
is seen that the function F(n*Q) is odd valued in e and, therefore, if 
the integral is to vanish for all n the function F(ris0) must be zero. 

The single equation specifying V(ti,Q) is then 
2 “ 

^ =f (ti) sin(e) , (H.17) 

dn 

where 

6^ = 2[l-cos(e)] = 4 sin^(e/2) . 

As before, the solution to the problem of vanishing stresses and 
displacements at infinity and uniform compression on the ends of broken 
fibers will now be sought. 

The boundary conditions are given by 

= 0 as ri “ 5 

for all fibers, 

, at n = 0 , (M.18) 

n = 0 , 

Equation (H.17) has the complete solution satisfying vanishing 
stresses and displacements remote from the damage as 

V(n,e) = A(e)e / sinh[6(n“t)]f^ (t)dt , (H.19) 

n 

where the functions A(e) and f-|(t) must be determined from the remaining 
boundary conditions. Using equation (H.14), the displacement is given by 


dV 

n _ — 

dn ^F 


= - 1 


for all broken fibers, and 
V^(0) = 0 at 


for all unbroken fibers. 
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Vn(n) =- / r A(e)e / sinh[6(n-t)]f, (t)dt 

" o L * n ' J 


n 

The remaining boundary conditions give 
dV„(0) 

- = t ^ { -i 5 A(e) - sin(e) . 

"" O 


sin(ne)d0 . 

(H.20) 


^U; o IT r ” T 

-3 - — f < - 6A(e) - sin{e) / cosh[6(n-t)]fn (t)dt > sin(ne)d0 = -1 , 

^ 0 I o • j 

(H.21) 

for all broken fibers, and 


V^(0) = ^ /^|^A{0) + sinh(6t)f-j (t)dt j sin(n0)d0 = 0 , 

° (H.22) 

for all unbroken fibers. Equation (H.22) is solved exactly by taking 


A(o) + / sinh(6t)f-| (t)dt = z B sin(me) , (H.23) 

0 ' m=l 

where are constants. Equation (H.22) then gives a system of 
N algebraic equations for N constants, B^ , in terms of f-|(n) which is 
as yet unknown. Eliminating A(e) between equations (H.20) and (H.23), 
the displacement Vj^(n) can be given as 

9 7T -(S ^ 

V(n)=— /G z B sin(me)sin(n0)d0 
n ^0 m=l 

- I fj f^(t) [c^(|t-n|) -Cj,(t+n) 

where 

c ( 5 ) = 2 / Sin(ne)d9 

n TT Q 0 

and 

6=2 sin(0/2) . 

Equation (H.21) then becomes 

2 IT f , / °° -^t , . 

— / ( -6 z B„ sin(m0) + sin(0) / e V-,(t)dt 
0 k m=l 0 

<» ^ 

- sin(0) / e T^(t)dt>sin(n0)d0 = -1, (H.25) 

0 ® j 

where f-|(t) has been replaced by T^-V-|(t) . 


dt , (H.24) 



Equations (H.IO) along with (H.24) gives 


V-| E sin(me)sin(e)d 0 


+ ^/°°{vi(t)-? 3 (t)jfc^{|t-n|)-C^(t+n)j dt , (H.26) 


where 




do 


Equations (H.25) and (H.26) must be solved simultaneously for V-|(ri) 
and the constants . Once Vi(n) and B^ are known, the displacement 
of any fiber n can be obtained from equation (H.20) and hence the fiber 
stresses. 

These equations were solved as follows: 


(I) An initial set of constants B^ was determined by setting f-j(t)=0, 
that is, V-|{t)= ^^(t) in equation (H.25) giving 


2 7T 

TT 0 


N 

-6 I 
m=l 


m 


sin(me)sin(n 0 )^ d 6 = - 1 


n = 1,2,...,N 


(H.27) 


(II) These initial constants were then substituted into the integral 
equation (H.26) and the function V-j(t) was determined for a known 
distribution of t,(t). 

a 

(III) Using V-j(t), a new set of constants B^ was computed from equation 
(H.25). 

(IV) This procedure was repeated until the unknowns changed less than 
a prescribed amount with the additional iterations. 

In the above solution, the unknown function, V-|(n) was assumed 
to be piecewise parabolic over the interval 0 5 n < of the form 


V’(n) = 


n + Yo 0 


i = 1 ,2, . . . ,k 
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when the interval was divided into k subdivisions. The function, V-j(n) 
contained 3k unknowns, (k+1) equations are obtained by requiring that 
the integral equation (H.26) be satisfied at the (k+1) end points, 
(2k-2) equations resulted from the requirement of continuity of the 
function V-|(n) and its first derivative between adjacent intervals 
and the last equation was given by setting V-j(n) = 0 as n Because 

the interval of integration was infinite, the piecewise linear approxi- 
mation for Vi(n) unlike in the preceding sections was found to be 
i nadequate. 


SOLUTIONS AND RESULTS 


Numerical solutions have been developed for all the specific prob- 
lems discussed in the previous sections, using an IBM 370 computer at 
Clemson University. For the case of a notch with no additional trans- 
verse or longitudinal damage, as considered in section A, the solution 
involves inverting the system of algebraic equations (A.21) for the 
constants, In sections B-G the inclusion of longitudinal damage 

requires the introduction of an unknown function g(n) representing the 
displacement difference between adjacent fibers in the damage zone. 
Satisfaction of the boundary conditions then gives a far more complex 
system of algebraic equations for the constants Bm, coupled with a 
Fredholm integral equation for the function g(n). The form of this pair 
of equations is identical for all sections B-G (see equations (B.46), 
(B.47), and (B.48)). The solution technique discussed in section B for 
these equations is appropriate for the remaining sections through G. 

For the half-plane considered in section H, even though no longitudinal 
damage is assumed, a set of equations similar to those of sections B-G 
still results. A discussion of an appropriate solution technique for 
the half-plane is given at the end of section H and some changes made 
necessary by unbounded intervals of integration are considered. 

For the present report the main emphasis will be to compare the 
results of the various models and demonstrate the agreement with 
available experimental results. More detailed results for the original 
Hedgepeth problem are given in [1]; the inclusion of longitudinal damage 
for one broken fiber is given in [2] and [3], and the general solution 
of section B is presented in [9]. 
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Some interesting observations concerning the case of broken fibers 
with no additional damage, as presented in section A and [1] are first 
considered, with particular emphasis on the behavior of the fiber stress 
in front of the notch. The normalized fiber stress in the first 
unbroken fiber was shown in [1] to have the following representation: 


^ F _ 4.6.8 (2r + 2) 

V a 375. 7 (2r + 1) 


( 1 ) 


where is the stress in the first unbroken fiber, and r is the total 

number of broken fibers. This equation may be written as 

,2r 

(2r + 1)] 


- ( 2 ) (r + 1) ! (r)! 


( 2 ) 


Using Stirling's formula for the asymptotic representation of (r)! when 
r is large, i.e. 

r 

( r ) ! /2Ir (^) , 


and substituting into equation (2), gives 

/n 


K = ^ /r 
r 2 


(3) 


If the total number of broken fibers is odd and the notch is symmetric 
about the y-axis, as is the case in all the full -plane problems consid- 
ered here, then r = 2N + 1, where N is the last broken fiber, so that 
equation (3) may be written as 


/n 




N+1 


+ 1 = y/| /N , for T 


2 ' V 2 ' ” > ' (4) 

The present failure criterion is similar to the "point stress" cri- 
terion presented by Nuismer and Whitney in £15], at least when applied 
to the models of sections A and B. That is, if transverse damage is not 
admitted then laminate failure is assumed to occur when the stress in 


the first unbroken fiber reaches the unnotched ultimate stress. The 



remote failure stress, a , is then given by equation (1) for any number 

oo 

of broken fibers, r, where is replaced by the unnotched strength 

r o 

Following Il5], a fracture toughness parameter can be defined as 
Kq = where a is the crack half-length. Kq is given by 

equation (4) for large crack lengths as 

Kq = > (5) 

where d is the fiber spacing, a - Nd, and /2/^N . 

A comparison can also be made between the solution of section A 
and that of a homogeneous isotropic plate with a central notch of length 
2 a {Griffith crack). The stress distribution at the crack tip along the 
x-axis is given by 



a 

OO 



X ^ a 


( 6 ) 


This equation can be modified for a unidirectional lamina in terms of 
the discrete fiber spacing and fiber index as 


nd 


" /(nd)^ - (N + - (N + C^j) 




, n N + 1, (7) 


where 



d 

n 


N 

Cn 


is the normalized stress in the fiber n, 
is the fiber spacing, 
is the fiber index 

is the index of the last broken fiber, and, 
is a constant needed to fit the data for each particular 
initial notch length specified by N. Note that (N + C|y])d 
represents an equivalent notch half-length. 
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For n = N + 1 this equation is set equal to equation (2), with r = 2N + 1» 
in order to determine the constant C|ij for any particular N. The 
following table of values for as a function of N results: 

N 0 1 10 1000 


Cfj 0.6614 0.6742 0.6805 0.6814 0.6815 



For large N equation (7), with n = N + 1, and equation (4) in place 
of equation (2) gives the limiting values of C|yj as 



from which 

Cm = 1 - - = 0.6815 = C. (8) 

IN TT 

It is then seen that Cfg is approximately equal to a constant, C, 
independent of N. Now, using the above values for C|yj the stress in 
fibers n > N + 1 in front of the notch was computed from equation (7) 
and compared with the solution of section A. For all notch lengths the 
agreement was found to be excellent as shown in Figure (13). Therefore 
the central notch in a unidirectional laminate using the shear-lag 
model has a square- root type stress distribution with an equivalent 
notch length given by ( N + Cf»j)d where the constant C|g - C = 1 - 1 /tt 
for all N. That is, the stress in the first unbroken fiber is given by 
the modified Griffith solution, equation (7), which also correctly 
predicts the stresses in the remaining unbroken fibers. 
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An analogous investigation concerning the maximum fiber stress and 
the fiber stress decay away from the notch when longitudinal and trans- 
verse damage is present, as given in sections B and C, is now discussed. 
The results are considerably different from the above and it is indi- 
cated that a basic assumption of [15] in which a square-root behavior 
was assumed to exist may be in error. That is, for cases having 
longitudinal and/or transverse damage the notch tip fiber stress as 
determined from sections B and C does not have a square-root form. To 
demonstrate this for longitudinal damage we used equation (7), with 
n = N + 1, to represent the fiber stress in front of the notch and 
determined the constant C|^ by comparing with the numerical results. A 
significantly different value of was found for each particular N as 
indicated in the following table. In constructing this table the ulti- 
mate stress in fiber N + 1 was taken as constant for all values of N; 
the remote stress and damage length, a, then changed with N. 


N 0 1 3 14 


Cfj 0.6459 0.5008 0.3437 -0.2364 


Further, using a specific value of N and the corresponding Cfj, the 
fiber stress computed for n > N + 1 from equation (7) was found to be 
different from the numerical results. An example of the fiber stress 
for seven broken fibers, and longitudinal damage consisting of yielding 
given by a = 5.0 is given in Figure (13). First, comparing with the 
results of the no damage, it is seen that the stress distribution is 
greatly reduced and distributed more uniformly to the adjacent fibers. 
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Second, the results are seen to be different from equation (7) but not 
drastically so. 

When transverse damage is included the stress field, shown in 
Figure (13) for four damaged fibers, is fundamentally different from 
that of the modified Griffith form, equation (7). A value of was 
not determined. It is indicated below that both longitudinal yielding 
and transverse damage are necessary in order to agree with the limited 
experimental results for boron/aluminum laminates. 

The failure criterion for longitudinal damage alone is the same as 
in section A. When the stress in the first unbroken fiber reaches the 
unnotched laminate stress and fails, then the stress in the adjacent 
unbroken fiber will be above the critical value. 

If transverse damage is also present the failure criterion is more 
complicated. The results from the model developed in section C indi- 
cated that stable transverse fiber breaks can occur. Failure of the 
first fiber does not then necessarily signify lamina failure as it may 
well require an increased applied stress to critically stress the next 
fiber. Recall that in the lamina, the fiber breaks occur in front of 
but not necessarily in line with the initial notch, with the matrix 
material remaining intact. Shear stress transfer from the matrix as 
well as some small axial matrix stress then gives a reduced load carry- 
ing capacity to the damage region in front of the notch as accounted for 
by the stiffness coefficient, y . Typical results are shown in Figure 
(14) for an initial notch length corresponding to seven broken fibers 
and a fiber stress in the damaged zone corresponding to a stiffness 
coefficient of y =0.90. It is seen that the transverse damage is 



stable until about seven additional fibers are broken. At this point 
the curve levels off and successive breaks require no increase in 
applied stress, An appropriate failure criterion is then to 

determine the applied stress, , giving unstable transverse damage. It 
is significant that the number of stable fiber breaks was approximately 
constant, independent of initial notch length. This constant damage 
zone size was about seven broken fibers which is in good agreement with 
the observed results of Awerbuch and Hahn [13]. 

The value of y , which represents the stress carried by the 
transverse damage region, was found to be that giving the best fit to 
the COD and strength curves along with a constant damage zone length. 

As the initial notch length increased, the load carried by the transverse 
damage region decreased, while the size of the region remained constant. 
Extensive experimental studies to investigate this behavior are 
underway. 

The results of the present methods are shown below to predict 
accurately the fracture behavior of both brittle and ductile matrix 
composites in terms of crack opening displacement (COD) and fracture 
strength. It then seems reasonable to conclude that the manner of 
stress distribution in the fibers is given with the same degree of 
accuracy. If this is true then it follows that, for the notched 
laminate with longitudinal and transverse damage, the fiber stresses in 
front of the notch are less severe than a square-root behavior. The 
failure criteria suggested by Nuismer and Whitney [15], either the 
"point stress" or the "average stress" criterion both assume a stress 
field in the front of the notch as 
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X 





(see equation (10) of[l5]). 


where a is the half notch length. Based on the above observations these 
criteria are perhaps not valid for unidirectional laminates. It 
appears then that longitudinal and/or transverse damage reduces the 
influence of a notch to a greater degree than an equivalent notch length 
model, with a square-root behavior, would predict. It is important to 
note that, even though the present study indicates the true fiber stress 
distribution, section C, to be different from [15j, results using the 
classical square-root behavior do predict the fracture characteristics 
of laminates with surprising accuracy. For example, see the work of Poe 
and Sova presented in references [16] and [17]. This is perhaps fortui- 
tous as the stress is not drastically different from a square-root form 
even though only the undamaged laminate of section A agrees accurately 
with the square-root behavior of equation (7). 

It has been demonstrated in [9] that for composites with a brittle 
matrix such as epoxy, the extensive splitting and the eventual insta- 
bility in the longitudinal direction can be predicted by introducing 
longitudinal damage alone as covered in section B. Results for cases of 
one and seven broken fibers are presented in Figure (15). 

For a composite with a ductile matrix such as aluminum some negli- 
gible splitting has been observed 114] and [1], but the stress in the 
fibers continues to increase with increased remote stress until the 
ultimate fiber stress is reached. In order to predict the fracture 
behavior of composites with a ductile matrix, such as boron/aluminum, it 
was found to be necessary to improve the model presented in section B. 



Effectiveness of these improved models, which are developed in sections 
C-E, will now be discussed. 

In all these cases, the model parameters such as number of con- 
strained fibers, stiffness coefficient, strain-hardening ratio, and 
constraint ratio were varied and their effect on the fracture charac- 
teristics of the laminate was investigated. A comparison of predicted 
values with the experimental results of [13] was done with respect to 
two fundamental characteristics, crack opening displacement and laminate 
strength. 

As explained in section B, the fiber axial displacement and 
stresses can be computed for different material properties, applied 
stress, and notch lengths. For a particular laminate a plot of COD vs 
applied load was obtained and then, by comparing with the experimental 
study of Awerbuch and Hahn in [13] for the unidirectional boron/ 
aluminum laminate, appropriate values of and the stiffness constant 
Gjvj/h were determined. The laminate used in [13] had the following 
material and geometric properties: 

Ef = 475 X 10^ Pa, 

Ap = 1.59 X 10-8m2 (o = 0.1427 mm), 

t = 0.165 mm/ply, eight plies 

^ul t = 3.98 X 10^ Pa 

w = width = 25.4 mm and 

fiber center line spacing = 0.178 mm. 

For a laminate having seven broken fibers, which corresponds to a 
notch length of about 1.27 mm, the load vs COD curve required specific 
values of and G^/h in order to give a "best fit" to the experimental 
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results of [13]. The numerical values of G)v|/h and determined in 
this manner are given in the following table for the various models: 

longitudinal transverse constraint strain 

damage, B damage, C layer, D hardening, E 

Gi^/h 65.4xlO^^N/m^ 115xlO^^N/m^ 148xlO^^N/m^ 35.9xlO^^N/m^ 

109xl0^N/m^ 88xlO^N/m^ 82xlO^N/m^ 164xlO^N/m^ 

For the strain-hardening case the modulus ratio was taken as Gr = 0.5 
and for the constraint layer Cr = 0.05. Using the above tabulated 
values in each model gave essentially the same COD curves for both seven 
and twenty-nine broken fibers with both being in dose agreement with 
the experimental results of [13]. This comparison is shown in Figure 
(16). 

The differences in the four models in predicting COD curves is then 
small although the particular values of the material constants Gf^/|/h 
and x^ are considerably different as seen from the table. 

Now using the above constants for G|v]/h and x^, the normalized 
ultimate fiber stress 



was computed for a specific ultimate fiber stress. Then the remote 
stress required to give a particular ultimate stress in the first 
unbroken fiber, or unstable extension in the case of section C, was 
determined from a plot of applied stress, a^, vs fiber stress, Op , Figures 
(17) and (18). This gives the strength of the laminate as a function of 
number of broken fibers (or notch length). This was repeated for all 
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the models in sections B through E. Figure (19) shows the plot of 
applied stress as a function of number of broken fibers for the differ- 
ent models along with the experimental results of [13], 

From Figure (19) it may be observed that even though COD was insen- 
sitive to model parameters, the strength curve was very much dependent on 
them. All the models predict the same trend in strength curves as those 
given by the experimental results. However, for the case of longitudi- 
nal damage, the model gives a much larger decrease in strength for small 
notch lengths than the experiment. At longer notches the model predicts 
the failure stress reasonably well. An increase in the fiber ultimate 
stress tends to rotate the curve about the knee of the curve without 
much change in the behavior for the shorter notches. 

Addition of a constraint layer does not make a significant differ- 
ence in the predicted strength curve. Variation in the constraint ratio 
had a similar effect as that of increasing the fiber ultimate stress in 
the longitudinal damage model. 

The strain-hardening model predicts strength values much lower than 
those of the longitudinal damage model, making it more and more notch 
sensitive. The best predicted values of the strain-hardening model are 
no better than those of the longitudinal damage model. 

The transverse damage model, however, predicts strength values very 
close to the experimental results and seems to account for an essential 
damage mode. By varying the transverse stiffness coefficient and number 
of constrained fibers, it was possible to move the predicted strength 
curve such that it fit the experimental results exceptionally well. One 
important observation was that the number of constrained fibers required 
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to fit the experimental data was approximately a constant at all notch 
lengths, indicating that the transverse damage zone size ahead of the 
notch tip was independent of the notch size. 

Introduction of transverse damage in addition to longitudinal 
damage then not only predicts COD in good agreement, at the same time it 
predicts strength values exceptionally well. The above computed value 
of is close to the yield stress for a homogeneous aluminum specimen. 
The equivalent stiffness, Gw/h, for h = 1.78 x 10“'^ m (i.e., the 
center-line distance between fibers) gives a shear modulus of 
G[/\ = 20,5 X 10^ |\|/m^ which is also close to the shear modulus of 
aluminum. For annealed aluminum the yield stress is about 95 x 10® N/m^ 
and the shear modulus is 27 x 10^ N/m^. 

Also of interest is the extent of the longitudinal yielding which 
is given by 

1 /EpAph 

L = = J Qi t ^ ^ fiber spacings. 

For the case of seven broken fibers the damage length corresponds to 
a = 11.5 or, therefore, L = 40.7 fiber spacings (approximately 6 times 
the crack length). For twenty-nine broken fibers a = 35, and L = 123.9 
fiber spacings (approximately 4 times the crack length). Experimental 
determination of the extent of longitudinal yielding is now being 
investigated by the writers; however, this damage does seem reasonable. 

With the above observations in mind it should be emphasized that in 
order to predict the fracture behavior of unidirectional composites with 
a ductile matrix it was essential to consider the transverse damage of 


matrix and fibers along with the longitudinal matrix damage as modeled 
in section C of this report. 

Note that the above results indicate an approximately constant 
transverse damage zone size of seven fibers, independent of initial 
notch length, which is in agreement with both Awerbuch and Hahn [13] and 
Nuismer and Whitney [15]. However, the disagreement with Nuismer and 
Whitney [15] as discussed above is that the present solutions predict 
fiber stresses in front of the notch to be different from and less 
severe than a square-root behavior. 

The final results to be presented in this report concern the solu- 
tions developed in sections F and G for the rectangular and circular 
cut-outs. Of principal interest will be the comparison of these results 
with the corresponding solutions for the notch. Extensive results are 
presented by Jones [18] where it is demonstrated that very little dif- 
ference exists between the three solutions concerning maximum fiber 
stress and the stress distribution in front of the damage. 

In Figure (20) results are presented for eleven broken fibers in 
the instance of no additional damage. The differences are largest in 
this case; any longitudinal or transverse damage brings the solutions 
closer together. Figure (21) gives the results for the same geometry 
with longitudinal damage, a = 4, and transverse damage of four fibers. 

In both cases the material constant, Hq, needed in the solution for 
the circular hole, was taken as 5.0 which corresponds to a volumetric 
ratio of 50% and Ep/G^ = 50. Within the range of realistic material 
properties for fiber dominated laminates, changes in Hq have little 
effect on the solution. 
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The behavior depicted by these two figures is typical of the gen- 
eral results and one can conclude that the shape of the initial damage 
is less significant in determining the laminate fracture properties than 
the total number of fibers broken. These results are fundamentally dif- 
ferent from the case of an isotropic plate in which the sharp notch has 
singular stresses near the notch tip with a square-root singularity; the 
rectangular hole also has singular stresses but with a power of less 
than one- half; and the circular hole has a stress concentration of 
three, independent of hole size. 

Future work in the area of this study will include a detailed 
experimental investigation of the damage growth in notched boron/ 
aluminum unidirectional laminates and the extension of the above half- 
plane solution to include damage and the presence of an adjacent half- 
plane of different material properties. 
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Figure 3. Two-dimensional unidirectional lamina with broken fibers 
and longitudinal matrix splitting and yielding (first 
quadrant) . 
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Figure 4. Two-dimensional unidirectional lamina with broken fibers, 
longitudinal matrix splitting and yielding and transverse 
fiber and matrix damage (first quadrant). 
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Figure 5. Boundary conditions along the x-axis for the transverse 
damage model . 
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Figure 6. Free-body diagram of an element consisting of the main 
lamina and the constraint layer. 
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Figure 7. Two-dimensional unidirectional lamina with a rectangular 
notch, longitudinal matrix splitting and yielding and 
transverse damage (first quadrant). 
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Figure 8. Two-dimensional unidirectional lamina with a circular 
cut out, longitudinal matrix splitting and yielding and 
transverse damage (first quadrant). 
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Figure 10. Assembly of boundary elements for the circular cut-out 
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Figure 12 


Free-body diagrams 
plane problem. 
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modified Griffith soln. and present soln. for no 
longitudinal or transverse damage 


present soln. for longitudinal damage alone 
(a = 5.0, S = O.O) 


modified Griffith soln. (equation (7)) for 
longitudinal damage alone (a = 5.0, B = 0.0) 


present soln. for longitudinal and transverse 
damages (a = 5.0, 4 constrained fibers) 
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Figure 13. Comparison of notch tip stress distribution for seven 
broken fibers. 









Figure 16. Comparison of present results (longitudinal damage model) 
with experimental study for a 25.4 mm by 8 ply laminate. 
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Figure 18. Maximum fiber stress in the first intact fiber at 
n = 0.0, as a function of applied stress for the 
case of no splitting and 29 broken fibers. 
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Number of broken fibers (2N+1) 
[crack length = fiber spacing x (2N+1)] 


Figure 19. Strength curve: Applied stress as a function of number of broken fibers. 
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